**Input the shapefiles: 

*Read in the shape-files and create the polygons and coordinates datasets: 

*Set the working directory: 
cd "$root/Data/Original"

*For the community borders: 

spshape2dta "gemeinde_grenze.shp", replace saving(comm)
*For the lake polygons: 
spshape2dta "k4seenyyyymmdd11_ch2007Poly.shp", replace saving(seen)


*Generate a new shapefile for the maps: 
spshape2dta "PLZO_PLZ.shp", replace saving(PLZ)

*Generate the label file for the city markers on the map: 
use PLZ.dta, clear 

bys PLZ: gen nb_obs=_n
gen PLZName=""
replace PLZName="Bern" if PLZ==3013
replace PLZName="Thun" if PLZ==3600
replace PLZName="Interlaken" if PLZ==3800 & nb_obs==2
replace PLZName="Biel" if PLZ==2502
replace PLZName="Burgdorf" if PLZ==3400
replace PLZName="Köniz" if PLZ==3098 & nb_obs==2

keep _ID _CY _CX PLZ PLZName
compress
save PLZ_labels, replace
